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INTRODUCTION 

Measures  of  stiffness  and  strength  of  most  structural  materials  are 
known  to  exhibit  considerable  spatial  statistical  variation.  The  inherent  dif¬ 
ference  in  the  behavior  of  tw  o  test  coupons  extracted  from  the  same  steel  bar 
is  but  a  typical  example.  These  differences  are  more  conspicuous  in  comoosite 
materials  such  as  structure,  ccncrjte  and  nsodcrn  fiber-reinforced  composites. 

The  random  nature  of  the  resistance  together  with  uncertainties  involved  in 
the  estimation  of  the  loading  conditions  made  it  necessary  for  the  structural 
engineer  to  adapt  the  so-called  "specified  minimum  strength"  that  could  lead 
to  an  uneconomical  design.  The  specified  minimum  strength  depends  only 
on  the  material  regardless  of  the  size  of  the  structural  member  to  be  designed. 
Accordingly,  the  variation  of  the  overall  resistance  with  the  size  of  a  structural 
member  cannot  be  predicted. 

The  weakest  link  hypothesis  was  introduced  to  explain  the  influence 
of  the  size  of  a  structural  member  on  the  magnitude  and  dispersion  of  strength.1-2 
Classically,  this  hypothesis  is  associated  with  "links"  having  statistically  indepen¬ 
dent  strengths;  this  overly  simplifies  reality  and  limits  the  scope  of  its  application. 
The  development  of  an  appropriate  criterion  for  the  variation  of  resistance  with 
size  for  prittle  materials  would  not  only  be  beneficial  in  structural  modeling,  but 
would  also  provide  the  means  for  a  better  interpretation  of  the  scatter  observed 
in  test  results  and,  thereby,  lead  to  better  design. 

The  overall  objective  of  this  study  is  to  establish  a  methodology  to 
deal  with  uncertainties  arising  from  the  randomness  of  the  properties  of  struc¬ 
tural  materials  with  particular  emphasis  on  structural  concrete.  To  this  end,  a 
new  and  more  general  interpretation  of  the  statistical  strength  distribution  is 
proposed;  this  interpretation  includes  the  classical  interpretation  as  a  special 
case.  To  illustrate  the  applicability  of  this  approach  to  composite  materials, 
a  probabilistic  model  was  constructed  for  the  digital  simulation  of  concrete 
strength  applicable  even  in  the  range  of  nonlinear  inelastic  response.  This 
model  is  compatible  with  the  finite  element  analysis  and,  therefore,  is  useful 
in  immediate  applications  such  as  in  planning  experiments,  in  parametric 
studies  involving  pertinent  strength  variate,  in  structural  modeling,  and 
in  analysis  and  design  of  structural  systems. 
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Statistical  Size  Effect 

Classical  Approach.  For  simplicity,  let  the  effect  of  statistical 
anisotropy  be  ignored  and  consider  the  statistical  scattei  of  a  measure  of 
strength,  denoted  by  X„,  of  a  hypothetical  reference  volume  element,  dv, 
of  a  microscopic  magnitude.  The  quantity  X0  may  be  interpreted  as  the 
strength  at  a  point  within  the  material. 

According  to  the  weakest  link  hypothesis,  if  Xn  denotes  the  strength 
of  an  aggregate  of  a  material  consisting  of  n  such  microscopic  volume  elements, 
then  X„  has  to  be  the  minimum  among  the  strengths  of  these  volume  elements 
constituting  the  aggregate  considered.  Under  the  assumption  that  these  strengths 
are  identically  and  independently  dislriDjicd  with  the  distribution  function 
Fv  (a),  the  distribution  function  F-  ?,d)  of  X„  is 

*o  "n  n 

Fx  (a)  =  1  -  fl  -  F  (a)]"  (1) 
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where,  by  definition,  the  distribution  function  Fx(a)  of  the  random  variable 
X  is  equal  to  the  probability  that  X  will  be  less  than  or  equal  to  a. 

As  n  increases,  ndv  =  v  becomes  finite,  and  Xv  {written  for  Xn) 
represents  the  strength  of  a  piece  of  material  of  a  finite  volume,  v.  It  is  usu¬ 
ally  postulated  that  the  distribution  function  of  X.,  approaches  tf  e  Weibull 
distribution,  which  is  one  of  the  asymptotic  distribution  functions  of  smallest 
values,  of  the  form 


FXy(a)  =  1  -  exp£-v(J)j  (2) 

where  A  and  K  ate  positive  constants  depending  only  upon  the  material 
properties. 

The  analytical  form  of  Equation  2  is  such  that  the  distribution 
function  of  the  strength  associated  with  volume  vt  is  located  to  the  left  of 
tha»  associated  with  volume  v2  if  v,  >  v2  as  schematically  shown  in  Figure  1 . 
This  indicates  the  statistical  size  effect:  an  aggregate  of  a  smaller  volume 
exhibits,  on  the  average,  a  larger  strength  than  that  of  a  larger  volume,  as 
long  as  the  two  aggregates  consist  oj  the  same  material. 

The  validity  of  such  a  postulate  depends,  among  other  things,  on  the 
mutual  independence  of  the  strengths  of  reference  volume  elements. 
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Strength  Argument,  a 


Figure  1.  Effect  of  size  on  strength  distribution  function. 


Proposed  Approach.  This  section  describes  the  proposed  interpreta¬ 
tion  of  the  statistical  strength  distribution  and  size  effect;  this  interpretation 
includes  the  classical  approach  as  a  special  case.  The  essence  of  the  proposed 
approach  lies  in  the  interpretation  that  the  strength  is  a  multidimensional 
random  process. 

For  simplicity,  consider  first  the  case  Of  a  one-dimensional  strength 
variation,  that  is,  the  variation  of  strength  along  a  line  arbitrarily  chosen  within 
the  material.  Without  loss  of  generality,  let  the  x-axis  represent  this  line  and 
XD(x)  denote  the  strength  of  the  material  at  a  point  foi  which  the  coordinate 
is  x.  The  proposed  approach  regards  XD(x)  as  a  one-dimensional  random  pro¬ 
cess.  One  may  visualize  X„(n)  as  the  axial  strength  variation  of  a  thin  rod  in 
which  the  radial  strength  variation  is  negligible. 

On  the  basis  of  experimental  evidence,  it  appears  reasonable  to  assume 
that  XD(x)  is  a  homogeneous  Gaussian  process  with  a  mean  value  Xm  and  a 
deviation  X(x)  from  Xm ;  hence,  XD(x)  can  be  expressed  as  follows; 

XJx)  =  Xm  +  X(x)  (3) 


Therefore,  the  standard  deviation,  denoted  by  o0,  of  X0(x)  equals  that  of  X(x). 

Consider  the  length  L  of  the  thin  rod  under  consideration  to  be  much 
larger  than  a  reference  length,  d.  Suppose  that  the  spectral  density  function 
o^g(k)  and  the  autocorrelation  function  o*  r(£)  of  the  process  X(x)  take  a 
normalized  form  such  as  that  shown  in  Figure  2;  that  is: 
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/hers  i  |n  •■$•••  :  ’imuni  ws1  e  number,  £  is  the  distance  between  two 
ooints  alu.  v '•  'r  jd,  and  Xm  {=  2i /\m)  is  the  wavelength  associated  with 
,<m.  The  •  ingthX„,  can  serve  as  a  correlation  distance  or  a  distance 
along  the  x-axis  over  which  appreciable  autocorrelation  occurs. 

According  to  the  proposed  approach,  the  ratio  between  the  correlation 
6  '  'ce  Am  and  the  reference  length  d  plays  a  major  role  in  the  interpretation 
cf  the  statistical  strength  variation  and  size  effect.  Figure  3  depicts  a  sample 
of  the  sire  jth  variation  XD(x)  with  X(x)  having  a  correlation  distance  con¬ 
siderably  shorter  than  the  reference  length  (Xm/d«1),  •  .nereas  another 
sample  of  X„(x)  with  X(x)  possessing  a  correlation  distance  considerably 
longer  than  d  (Am/d»1)  is  schematically  shown  in  Figure  4.  Asa  limiting 
case,  a  situation  can  be  considered  where  the  correlation  distance  approaches 
i.sro  and  the  autocorrelation  takes  the  form  6(£),  where  6(£)  is  the  Dirac 
delta  function.  This  corresponds  to  the  classical  approach  of  independent 
strength  variation.  Another  limiting  case  produces  a  complete  coherence  in 
which  the  correlation  distance  is  infinite.  Consequently,  no  strength  variation 
is  observed  along  the  x-axis  within  one  specimen;  however,  the  strength  still 
may  vary  from  one  specimen  to  another. 

Aside  from  the  two  limiting  cases  mentioned  above,  the  proposed 
approach  permits  correlation  among  the  strengths  of  various  elements.  The 
trend  of  such  strength  correlation  is,  in  general,  at  least  similar  to  that  of  the 
autocorrelation  function  of  X(x).  In  the  thin  rod  example,  if  the  correlation 
distance  Am  is  much  shorter  than  the  reference  length  d,  no  significant  corre¬ 
lation  is  expected  between  even  the  neighboring  elements;  whereas,  if  Am  is 
much  longer  than  d,  significant  correlation  is  expected  over  a  number  of  suc¬ 
cessive  elements  whose  total  length  is  approximately  equal  to  Xm .  Accordingly, 
for  the  purpose  of  simulating  the  strength  variation,  the  selection  of  the  param¬ 
eters  defining  the  process  X0(x)  may  be  guided  by  the  fact  that  the  spectral 
density  and  the  autocorrelation  functions  form  a  Fourier  transform  pair. 

Hence,  the  spectral  density  is  flatter  or  wide-banded  (more  concentrated 
or  narrow-banded)  if  the  correlation  distance  is  shorter  (longer). 
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Figure  2.  Graphs  of  spectral  density  and  autocorrelation  functions 
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To  interpret  the  statistical  size  effect,  consider  again  the  example  of 
a  thin  rod  for  which  the  strength  variation  in  the  radial  direction  is  negligible 
when  compared  with  that  in  the  axial  direction.  Invoking  the  weakest  link 
hypothesis  together  with  the  proposed  interpretation  of  the  strength,  the 
strength  of  a  rod-element  of  a  length  d  is  then  the  minimum  value  of  X„(x) 
or  Xm  plus  the  minimum  value  of  X(x)  in  that  element.  Therefore,  the 
strength  of  a  larger  element  is  likely  to  be  smaller  when  the  size  effect  is 
taken  into  account. 

The  preceding  idea  of  the  statistical  size  effect  can  easily  be  expanded 
into  more  realistic  two-  or  three-dimensional  situations.  Consider,  for  example, 
the  simple  beam  shown  in  Figure  5  where  the  thickness  c  is  small  enough  that 
the  statistical  variation  of  the  strength  across  the  thickness  is  negligible.  If  the 
strength  at  any  point  (x,y)  is  treated  as  a  function  X0(x.y),  then,  according 
to  the  weakest  link  hypothesis,  the  strength  of  an  element  of  the  beam,  such, 
aseof  area  Ae  =  dd'  (Figure  5),  is  the  minimum  value  of  XD(x,y)  within  the 
area  Ae. 


Thickness  -  c 


Figure  5.  A  simple  beamrgrid  for  strength  variation. 


It  again  appears  reasonable  to  assume  that  X0(x,y)  is  a  homogeneous 
(although  two-dimensional  this  time)  Gaussian  process  consisting  of  the  n.^an 
value  Xm  and  the  deviation  X(x,y)  from  Xm  : 

X0(x,  y)  =  Xm  +  X(x,y)  (5) 
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where  X(x,  y)  is  a  homogeneo  us,  two-dimensional  Gaussian  process  with 
mean  zero,  standard  deviation  a0.  and  "normalized"  autocorrelation  func¬ 
tion  !•(£,, $2 ).  The  variables  £,  and  £2  indicate  the  separation  between  two 
points  (x,,  y, )  and  (x2,  y2):  £1  =  x2  -  x,  and  £2  =  V2  "  Vi  Since  X(x,  y) 
is  homogeneous,  the  double  Fourier  transform  q(kv  k2)  of  r(£1#£2)  of  the 
form 


OO  po 

9<kvM  =  /  /  r^i^2>eH<kltl  +  k2t2,d£,d£2 


is  real  and  n  jnnegaiive  by  virtue  of  Bochner's  theorem.3 

The  ixpansion  of  the  same  idea  into  the  three-dimensionai  case 
involving  Xc(x,  y,  z)  is  straightforward. 

Digital  Simulation  of  Strength  Variation 

The  strength  variation,  which  is  interpreted  herein  as  a  random  process, 
is  simulated  utilizing  the  method  developed  in  Reference  4. 

Consider  first  the  one-dimensional  strength  variation  X0ix)  -  Xm  + 
X(x).  Since  Xm  is  a  constant,  X(x)  with  mean  zero  and  autocorrelation  o2  r(£) 
is  the  process  to  be  digitally  simulated. 

According  to  Reference  4.  if  G{x)  denotes  a  homogeneous  Gaussian 
process  with  mean  zero,  unit  standard  deviation,  and  autocorrelation  r(£),  then 
G(x)  can  be  t  xpressed  as 

GW  ■  (s) 

where  k;  (andom  variables  identically  and  independently  distributed 
with  density  function  g(k),  which  is  the  Fourier  transform 
of  r(£) 

M  =  .1  large  positive  integer  (say,  100) 

4>j  =  1  nndom  variables  identically  and  independently  distributed 
with  uniform  density  between  zero  and  2w  and  is  independent 
of  kj  (j  =  1,  2, . . .  M) 

Equation  7  can  be  utilized  for  the  digital  simulation  of  G(x)  through  the 
replacement  of  the  random  variables  kj  and  4>;  by  their  realizations.  Since 
X(x)  can  be  written  as 


£  cos(kjX  +  4>j)  (7) 
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X<x)  =  o0G(x) 


(8) 


It  is  evident  that  a  realization  of  the  minimum  value  of  XU)  in  an  interval  of 
length  d  is  equal  to  that  of  G(x),  in  the  same  interval,  multiplied  by  oQ. 

A  similar  technique  can  be  applied  to  simulate  the  strength  variation 
over  a  two-dimensional  domain.  In  this  case,  the  strength  is  interpreted  as 
the  process  X0U,y)  which  can  be  decomposed  into  the  mean  value  Xm  and 
the  deviation  XU)  as  given  by  Equation  5.  Following  the  same  line  of  reason¬ 
ing  as  in  the  one-dimensional  case,  a  two-dimensional  homogeneous  Gaussian 
process  GU.y)  can  be  simulated  Digitally  in  the  form 

/o\l/2  M 

GU.y)  =  £  cos(k1i  x  +  k2iy  +  4»j)  (9) 

where  kti  and  k2i  =  jointly  distributed  random  variables  with  identical 
joint  density  function  gik^kj)  irrespective  of  i; 
independent  of  k,j  and  k2j  (i  ^  j) 

M  =  a  large  positive  integer 

4»j  =  random  variable  uniformly  distributed  between  zero 

and  2 it;  independent  of  (i  *  j)  and  of  any  set  of 

k1}  and  k2j 

Replacing  k^,  k2i.  and  in  Equation  9  by  their  respective  realizations 
generates  a  realization  of  GU.y).  A  realization  of  X(x,y)  can  then  be 
determined  from  the  relation 

X(x,y)  =  o„  G(x,y)  00) 

Obviously,  realizations  of  X0(x,y)  can  be  obtained  from  those  of  X(x,y) 
utilizing  Equation  5. 

Again,  the  expansion  of  the  method  of  simulating  material  stiength 
into  the  three-dimensional  case  is  straightforward. 

Estimation  of  Parameters  From  Observed  Data 

This  section  deals  with  the  procedure  for  establishing  the  link  between 
the  results  from  the  digital  simulation  and  the  data  observed  in  laboratory 
experiments.  Such  a  link  would  provide  the  parameters  for  defining  the  ran¬ 
dom  process.  This  is  done  by  fitting  the  distribution  function  resulting  from 
the  assumed  random  process  to  that  constructed  from  the  observed  data. 
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Let  Fd(a)  denote  the  distribution  function  of  the  minimum  value  of 
the  process  XD(x)  in  an  element  of  length  d  (Figures  3  or  4).  Then,  Fd(a) 
depends  in  general  not  only  on  d  but  also  on  the  parameters  defining  the 
process  X0(x).  When  Xc(x)  is  a  homogeneous  Gaussian  process  with  mean 
Xfn ,  standard  deviation  a0.  and  autocorrelation  function  a\  r(£)  as  assumed 
in  this  study,  the  distribution  function  Fd(a)  depends  on  Xm .  0o ,  and  the 
analytical  form  of  r{£)  as  well  as  d.  The  procedure  for  establishing  such  a 
dependence  is  described  below. 

For  the  thin  rod  example,  select  an  autocorrelation  function  r(£)  for 
an  interval  of  length  L  of  the  rod  consisting  of  N  subintervals  d.  In  this  case, 

L  =  N  d  is  chosen  in  such  a  manner  that  N  is  a  reasonably  large  positive  integer 
(say,  200).  Simulate  the  corresponding  homogeneous  Gaussian  process  G(x) 
as  outlined  in  the  previous  section.4  The  minimum  value  of  the  realization  of 
G(x)  in  each  of  the  subintervals  d  constitutes  an  element  of  a  sample  of  size  N. 
Because  N  is  large,  the  distribution  function,  denoted  by  Fd(a),  of  the  mini¬ 
mum  value  of  G(x)  over  the  interval  L  can  be  inferred  from  such  a  sample 
with  good  engineering  confidence. 

As  stated  earlier,  the  strength  variation  X(x)  can  be  simulated  as 
oQ  G(x).  and  that  a  realization  of  the  minimum  value  of  X(x)  in  an  interval 
of  length  d  is  equal  to  that  of  G(x)  in  the  same  interval  multiplied  by  oD.  This 
indicates  that  the  corresponding  distribution  function  of  the  minimum  for 
XQ(x)  =  Xm  +  X(x)  can  simply  be  obtained  from  the  distribution  function 
F°(a)  of  the  minimum  for  G(x)  as  follows: 


■-m 


Obviously,  a  realization  of  G(x)  generated  over  the  interval  L  can  be 
utilized  to  estimate  ihe  distribution  function  Fd(a)  with  different  values  of 
d  as  long  as  L/d»  1 .  Accordingly,  in  applying  this  Monte  Carlo  method  for 
the  estimation  of  Fd(a),  only  one  realization  of  G(x)  of  sufficient  length  has 
to  be  generated  for  a  class  of  processes  with  the  normalized  autocorrelation 
function  identical  to  that  of  G(x). 

One  should  recall  that  the  distribution  function  Fd(a)  given  by  Equa¬ 
tion  1 1  was  derived  from  theoretical  'onsideralion  rather  than  experimental 
observa  ions,  and  that  herein  Fd(a)  r  spends  on  a  tentative  set  of  parameters 
lr(£),  Xm,a0,  and  d] .  7 he  select ion  of  the  appropriate  values  for  the  elements 
of  a  set  calls  for  a  laboratory  experiment  such  as  the  one  outlined  below. 

Consider,  again,  the  thin  rod  example  where  the  strength  is  measured 
for  N  subintervals  d  of  an  interval  of  length  L.  This  time,  however,  the  experi¬ 
ment  is  conducted  in  the  laboratory  for  an  "actual"  thin  rod.  This  experiment 


provides  a  sample  of  N  strength  values  taken  from  the  population  of  strength 
value?.  Fd(a).  From  the  results  of  such  an  expei  iment,  it  is  possible  to  estimate 
approximately  the  analytical  form  of  r(£).  This  makes  it  possible  to  find  Fd(a) 
by  the  method  of  digital  simulation  described  earlier.  The  remaining  two 
parameters  (the  mean  Xm  and  the  standard  deviation  oc)  can  then  be  esti¬ 
mated  in  such  a  way  that 

/a  -  Xm\ 

-  F‘,a)  <12> 

As  far  as  the  problem  in  hand  is  concerned,  the  validity  of  Equation  12  insures 
conformity  of  concept  to  reality. 

The  preceding  approach  can  readily  be  extended  into  a  two-dimensional 
case.  Recalling  the  example  of  the  simple  beam.  Figure  5,  the  distribution  func¬ 
tion  FA(a)  of  the  resistin';  strength  of  the  element  e  of  area  Ae  =  d  d*  can  be 
obtained  from  the  distribution  function  FA(a}: 


'.••here  F£(a)  is  the  distribution  function  of  the  minimum  value  of  G(x,y) 
within  the  element  e  of  dimensions  dd'c.  The  distribution  function  F£(a) 
is  to  be  inferred  on  the  basis  of  digital  simulation  of  G(x,  y)  over  a  two- 
dimensional  region  (for  example,  the  rectangular  area  Lh  of  Figure  5) 
consisting  of  a  large  number  of  subregions  (elements  d  d'.  Figure  5).  The 
autocorrelation  function  r(|1f  $2)  of  G(x,y)  is  required  for  such  a  simulation. 
Similar  to  the  one-dimensional  case,  the  parameters  Xm.  and  o„  can 

be  found  from  observed  results  of  strength  measurements  performed,  for  exam¬ 
ple,  on  each  of  N  elements  of  size  d  d'  c  of  a  rectangular  region  of  thickness  c. 

Again,  the  expansion  of  this  procedure  into  the  three-dimensional  case 
is  straightforward. 

With  an  additional  effort  of  performing  simple  transformations,  one 
can  treat  without  much  difficulty  the  strength  process  which  is  not  Gaussian. 
However,  in  the  face  of  the  limited  amount  of  information  available  at  the 
present  time,  such  an  effort  may  not  be  worth  the  expense. 

Finite-Element  Failure  Analysis 

This  section  demonstrates  that  it  is  possible  to  perform  a  more  realistic 
structural  failure  analysis  when  the  proposed  statistical  strength  model  is  utilized. 
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When  a  failure  analysis  is  sought,  a  failure  criterion  is  always  required. 
Failure  criteria  normally  call  for  comparison  between  a  computed  or  measured 
load-dependent  action  (for  example,  stress,  strain. . . .  etc.)  and  a  known  or 
assumed  ultimate  value  of  such  action.  This  ultimate  value  is  usually  referred 
to  as  strength  or  resistance.  In  what  follows,  it  is  assumed  that  (1 )  the  'oad- 
dependent  action  is  computed  by  means  of  the  finite  element  method,5  and 
(2)  the  resistance  is  determinable  from  a  process  such  as  XQ(x,  y)  for  which 
the  parameters  Xm,  oQ,  and  r(£1#£2)  are  estimated  (based  on  strength  mea¬ 
surements)  with  reasonable  engineering  accuracy.  Without  loss  of  generality, 
it  is  also  assumed  that  the  failure  criterion  is  applicable  for  actions  within  the 
finite  elements  rather  than  at  the  nodal  points  of  the  idealized  structure.  The 
ensuing  discussion  addresses  itself  mainly  to  the  question  of  determining  the 
resistance  of  each  of  the  finite  elements. 

Consider  a  two-dimensional  failure  analysis  of  an  arbitrary  structure. 
Figure  6.  for  which  the  strength  variation  is  negligible  in  the  direction  normal 
to  the  plane  of  interest.  For  the  purpose  of  discussion,  a  coarse  grid  of  quadri¬ 
lateral  finite  elements  is  shown  in  Figure  6;  the  selection  of  the  "best"  finite 
element  grid  is  considered  outside  the  scope  of  this  study. 


Figure  6.  Finite  element  idealization  of  a  two-dimensional  structure. 
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Utilizing  the  approach  outlined  earlier,  a  number  of  realizations  of 
X0(x,  y)  can  be  generated  digitally.  Each  realization  will  be  equivalent  to  a 
structure  chosen  arbitrarily  from  the  population  of  structures  with  the  specified 
form  of  statistical  strength  variation.  Obviously,  the  value  of  X0{x,y)  can  be 
evaluated  for  any  pair  of  finite  values  of  x  and  y.  However,  when  x  and  y  are 
restricted  to  the  region  representing  the  structure,  the  corresponding  values  of 
X0(x,  y)  become  of  special  significance  for  they  represent  the  spatial  strength 
variation  of  the  structure.  For  such  digital  simulation  of  strength,  an  arbitrary 
grid  may  be  selected,  and  the  value  of  Xe(x,  y)  may  be  determined  at  the 
intersection  points  of  the  grid.  In  what  follows,  the  grid  for  digitally  evalu¬ 
ating  values  of  X0(x,  y)  will  be  referred  to  as  the  strength  grid. 

For  the  special  case  in  which  the  strength  grid  coincides  with  that 
ct  istituting  the  finite  elements,  the  intersection  points  are  the  nodal  points 
of  the  finite  element  grid.  Thus,  the  number  ne  of  values  of  X0(x,  y)  for  points 
(x,  y)  falling  within  each  finite  element  is  equal  to  the  number  of  exterior  nodal 
points  for  that  element.  Applying  the  weakest  link  hypothesis  at  the  finite 
element  level,  the  strength  of  any  element  is,  then,  the  minimum  of  n,  values 
of  XD(x.y)  for  that  element.  The  strength  of  each  finite  element  thus  obtained 
can  be  compared  with  the  load-dependent  action  obtained  from  the  finite  ele¬ 
ment  solution  to  determine  whether  failure  at  the  elemr  nt  level  has  occurred. 

In  the  preceding  discussion,  the  strength  grid  was  assumed  to  coincide 
with  the  finite  element  grid.  However,  this  is  not  necessary  because  the  grid 
of  finite  elements  depends  mainly  on  the  expected  stress  gradient  within  the 
structure  under  consideration:  steeper  gradients  call  for  finer  finite  element 
grids.  Similarly,  the  selection  of  a  strength  grid  should  depend  on  the  mariner 
in  which  the  strength  varies.  More  precisely,  the  size  of  the  strength  grid  should 
be  in  harmony  with  the  autocorrelation  function:  shorter  correlation  distances 
should  call  for  finer  strength  gnds.  For  expediency,  however,  it  is  not  advisable 
to  make  the  strength  grid  coarser  than  the  finite  element  grid  regardless  of  the 
length  of  the  correlation  distance.*  The  reason  is  that  difficulty  will  arise  in 
defining  the  strength  of  any  finite  element  for  which  ne  is  zero. 

In  summary,  it  is  recommended  that  each  finite  element  contain  a 
strength  grid  such  that  the  nodal  points  of  that  element  also  will  be  the  inter¬ 
section  points  of  that  grid.  The  size  of  the  strength  grid  within  each  element 
should  be  made  finer  as  the  correction  distance  gets  shorter.  In  this  context, 
finer  grid  means  that  the  number  of  the  intersection  points  of  the  strength  grid 
exceeds  the  number  of  nodal  points  of  the  finite  element  as  shown  in  Figure  7. 


•  The  only  exception  is  the  limiting  case  where  the  correlation  distance  is  infinitely  long; 
in  this  case,  all  the  elements  have  the  same  strength,  and,  therefore,  the  construction  of 
a  strength  grid  is  irrelevant. 
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nodal  point 


As  a  typical  example,  consider  a  failure  analysis  of  a  two-dimensional 
concrete  structure.  The  strength  of  an  element  of  plain  concrete  of  a  certain 
size  under  a  two-dimensional  state  of  stress  can  be  expressed  in  terms  of  a  fail¬ 
ure  envelope.  Figure  8,  in  the  ot  -  o2  plane  of  principal  stresses.6  The  domain 
.n  the  a,  -  a2  plane  enclosed  by  such  an  envelope  represents  the  various  states 
of  stresses  that  can  be  imposed  without  causing  a  failure  of  that  element.  The 
equation  of  this  envelope  can  be  written  in  the  form 

R(x,y;0)  =  Xo(x,y)Ro(0)  {14} 

where  0  is  a  measure  of  the  stress  path,  Ro(0)  represents  a  reference  strength, 
and  R(x,y;0)  is  a  random  process  describing  the  statistical  variation  of  the 
strength.  Equation  14  implies  the  hypothesis  that  the  shape  of  the  failure 
envelope  is  geometrically  similar  but  its  size  varies  statistically;  this  hypothe¬ 
sis.  however,  is  not  essential  to  the  present  analysis.  According  to  this 
hypothesis,  the  relative  positions  of  the  failure  envelopes  for  two  intersection 
points P1(x1,y1)  and  P2(x2,y2),  Figure  7.  with  Xc(x1,y1)>  X0(x2,y2)  is 
shown  schematically  in  Figure  8.  If,  for  example,  the  reference  envelope  is 
associated  with  the  median,  the  failure  envelope  for  any  arbitrary  point  P{x.  y) 
has  a  probability  of  0.5  to  be  located  outside  (or  inside)  the  reference  envelope. 
This  also  implies  that  the  probability  of  X0(x,y)  for  any  point  P(x,y)  being 
larger  (or  smaller)  than  unity  is  0.5. 
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Failure  Envelopes 

. .  Reference  :  R(*.  y;  8)  =  R0(8) 

- P,(X1.  Vi»:  RUi.Vl:8)  =  X0<xv  v,)  R0{8> 

—  —  —  P2<*2-  V2,:  R<*2-  y2- 0>  “  xo(x2-  v2*  Ho19* 


Figure  8.  Envelopes  for  failure  of  concrete  under  biaxial  stresses. 


Assuming  that  the  beam  shown  in  Figure  5  is  made  of  plain  concrete, 
the  present  model  of  statistical  strength  variation  makes  it  possible  to  simulate 
X0(x.y)  and.  therefoie,  to  assign  spatially  correlated  failure  envelopes  to  the 
finite  elements.  As  the  central  load  is  increased  from  zero,  the  principal 
stresses,  from  the  finite  element  solution,  are  examined  until  "failure"  of  an 
element  :s  identified.  The  value  of  the  central  load  that  causes  this  condition 
is  recorded,  and  the  procedure  is  repeated  for  a  reasonably  large  number  of 
(simulated)  beams.  Accordingly,  one  can  estimate  the  statistical  characteristics 
of  the  structural  failure.  Obviously,  such  a  failure  analysis  is  extremely  useful 
in  the  prediction  of  structural  integrity.  It  is  important  to  note  that  the  stress 
analysis  has  to  be  performed  only  once  if  the  material  properties  such  as  the 
elastic  constants  are  not  statistical. 

In  this  example  of  a  plain  concrete  beam,  the  weakest  link  hypothesis 
was  assumed  to  govern  explicitly  at  the  finite  element  level  and  implicitly  at 
the  global  or  the  entire  structure  level.  This  implic't  assumption  at  the  giobal 
level  may  be  realistic  for  a  plain  concrete  beam  (brittle  failure)  but  not,  in 
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general,  for  a  reinforced  concrete  beam  (ductile  failure).  It  is  to  be  emphasized 
that  the  proposed  aDproach  suggests  the  weakest  link  hypothesis  only  at  the 
local  (finite  element)  level  but  does  not  require  such  assumption  at  the  global 
level. 

In  the  preceding  sections,  the  random  process  X0(x,y)  was  defined 
as  a  measure  of  the  strength  interpreted  as  the  resistance  associated  with  fail¬ 
ure.  It  is  equally  likely,  and  perhaps  as  important,  to  utilize  the  same  procedure 
of  digital  simulation  when  the  process  X0(x,  y)  is  defined  as  a  measure  for  any 
material  property  (for  example,  mouulus  cf  elasticity)  of  a  finite  element.  In 
such  a  case,  a  decision  has  to  be  made  as  to  whether  the  weakest  link  assump¬ 
tion  is  realistic  even  at  the  finite  element  level;  an  averaging  procedure  may 
provide  a  better  alternative. 

APPLICATION:  DIGITAL  SIMULATION  OF  CONCRETE  MODELS 

This  section  presents  a  numerical  example  of  applying  ‘he  proposed 
approach  to  investigate  the  effect  of  size  on  the  strength  variation  of  speci¬ 
mens  made  of  concrete.  Because  of  the  lack  of  experimental  data  concerning 
the  spatial  strength  variation  of  concrete,  assumptions  had  to  be  made  in  order 
to  define  the  strength  as  a  random  process.  However,  instead  of  being  abso¬ 
lutely  arbitrary,  such  assumptions  were  made  to  reflect  some  experimental 
evidence.  Also  discussed  is  the  influence  of  the  choice  of  the  strength  grid 
on  the  resulting  strength  of  the  simulation  models. 

Failure  Level  and  Simulation  Parameters 

Failure  of  a  structural  material  is  normally  associated  with  a  certain 
level  or  more  precisely  a  certain  dimension  at  which  failure  is  observed  or  esti¬ 
mated  (for  example,  microscopic  failure,  macroscopic  failure,  etc.).  It  is  general 
practice  to  consider  that  failure  of  a  structural  material  occurs  when  some  mea¬ 
sure  of  stress  (or  strain)  has  reached  or  exceeded  its  ultimate  value  (usually 
referred  to  as  the  strength).  By  definition,  stresses  are  associated  with  lengths 
or  dimensions.  But,  strains  are  dimensionless;  they  relate  displacements  to 
some  physical  dimension  or  gage  length.  However,  it  is  weil  known  that  the 
value  of  the  strain  is,  in  general,  a  function  of  the  gage  length.  In  other  words, 
strains  are  dependent  on  dimensions  in  spite  of  their  apparent  nondimensionality. 
This  paradox  normally  disappears  when  the  gage  length  is  standardized  whether 
this  is  done  to  relate  strains  either  to  stresses  or  to  an  ultimate  value  for  failure 
analysis.  Accord’ngly.  failure  is  associated  with  dimensions  implicitly  (or  explic¬ 
itly)  when  the  failure  criterion  is  based  on  strains  (or  stresses).  The  selection  of 
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the  dimension  with  which  failure  is  associated  is  limited  by  the  precision  of 
the  instruments  available  and  governed  by  the  consequences  of  failure  (for 
example,  collapse)  at  the  selected  dimension. 

In  considering  structural  concrete,  design  engineers  are  generally 
interested  in  avoiding  major  failures,  that  is,  the  failures  associated  with  large 
dimensions.7  In  this  case,  the  assumption  that  plain  concrete  behaves  as  a 
homogeneous  material  is  reasonable.  However,  this  assumption  becomes 
unrealistic,  and  composite  action  needs  to  be  accounted  for  when  failure  is 
associated  with  smaller  dimensions  as  in  the  case  of  bond  failure  between 
cement  paste  and  aggregate.8,9  In  considering  failure  at  a  much  smaller 
dimension,  propagation  of  concrete  cracking  is  thought  to  initiate  from 
capillary  pores  (voids  in  the  cement  paste  that  were  originally  filled  by 
water  but  have  not  become  filled  by  the  products  of  hydration).10 

Generally,  as  the  dimension  with  which  failure  is  associated  decreases, 
the  strengths  and  maximum  stresses  increase,  and  the  maximum  and  minimum 
values  occur  at  shorter  intervals.10  The  spatial  variation  of  a  measure  of  one- 
dimensionai  strength  of  concrete  as  influenced  by  the  gage  length  may  take 
the  form  shown  in  Figure  9. 

In  order  to  make  the  numerical  example  of  practical  value  to  designers, 
the  failure  considered  is  that  associated  with  dimensions  comparable  to  the 
size  of  standard  concrete  cylinders.  For  such  sizes,  the  most  important  param¬ 
eters  that  influence  the  variation  of  strength  with  size  are  the  maximum 
aggregate  size  and  the  fineness  modulus.11"13  Based  on  the  theory  of  dimen¬ 
sional  analysis.  Brown  proposed  consideration  of  the  maximum  aggregate  size 
in  studies  of  size  effect.14  Obviously,  many  parameters  such  as  the  aggregate 
gradation,  mix  proportions,  method  of  mixing,  placing,  curing,  age.  and  rate 
of  loading  influence  the  strength  of  a  concrete  specimen.  The  question  is 
whether  such  influences  change  with  the  sizes  considered. 

For  the  numerical  example  presented  below  the  following  assumptions 
were  made: 

1 .  Strength  in  any  direction  is  evaluated  from  failure  of  specimens 
with  a  minimum  dimension  of  A  such  that  A  =  4  Am  where  Am  is  the  maxi¬ 
mum  aggregate  size.15  Such  a  specimen  will  be  referred  to  as  the  reference 
specimen. 

2.  The  mean  strength  Xm ,  the  standard  deviation  o0.  and  the 
autocorrelation  function  r(£)  are  evaluated,  in  the  manner  discussed  earlier, 
from  specimens  conforming  with  Assumption  1  above: 

A  =  d  =  4Am 

and,  therefore, 

FJa)  =  Fd(a) 
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3.  The  autocorrelation  function  takes  a  form  similar  to  that  shown 
in  Figure  2  such  that 


\n  =  24  Am 


Description  of  Models 

For  a  concrete  miy  with  a  maximum  aggregate  size  Am  =1/4  inch, 
consider  a  reference  specimen  of  dimensions  A  Ac  inches  such  that  A  =  4  Am  = 
1 .0  inch,  and  assume  that  the  strength  does  not  vary  with  he  thickness  c.  The 
strength  of  the  reference  specimen  is  determined  by  applying  uniform  pressure 
along  two  opposite  edges  without  constraining  the  movement  of  any  edge.  Fig¬ 
ure  10.  The  following  parameters  are  assumed  to  have  been  determined  from 
testing  a  large  sample  of  reference  specimens: 

Xm  =  10,000  psi  a0  =  2,000  psi  Xm  =  6  inches 


The  proposed  approach  was  utilized  to  determine  the  variation  of 
strength  when  the  specimen's  dimensions,  in  inches,  vary  from  1  x  1  x  c  to 
3x  3xcand  to  12  x  12xc  for  the  same  type  of  loading.  Figure  10.  For  this 

loading  condition,  the  stress  distri¬ 
bution  does  not  vary  with  the  size 
of  the  specimen.  Accordingly,  the 
finite  element  grid.  Figure  11a,  con¬ 
sisted  of  the  same  number  of  finite 
elements  for  both  3  x  3-inch  and 
12  x  12-inch  specimens.  The  size 
de  x  u,  of  each  finite  element  was 

I  x  1  inch  and  4x4  inches  for  the 
smaller  and  the  larger  specimens, 
respectively.  The  strength  grid 
within  each  finite  element.  Figure 

I I  b,  was  selected  such  that  the  dis¬ 
tance  ds  between  the  intersection 
points  is  small  relative  to  the  corre- 

FiguretO.  Reference  concrete  specimen.  lat ion  distance  Xm.  This  was  done 

so  that  the  influence  of  the  fineness 
of  the  strength  grid  on  the  calculated 

strength  becomes  insignificant.  The  strength  distribution  function  for  the 
3  x  3-inch  specimens  was  evaluated  for  two  ratios  of  Xm/d,,  namely  6  and  24. 
When  the  distribution  functions  for  the  two  ratios  were  found  to  be  almost 


p  =  uniform  pressure 
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Am  »  1/4  inch 
6  inches 
Xm  -  10.000  psi 
oQ  »  2,000  psi 
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identical,  the  Xm/dt  ratio  of  6  was  selected  for  the  sttergth  evaluation  of  the 
12  x  12-inch  specimens.  Table  1  lists  the  properties  of  the  specimens  for  which 
the  strength  distribution  functions  were  evaluated. 


(a)  Finite  element  grid  for  each  specimen. 


(b)  Strength  grid  for  each  finite  element. 

Figure  11.  Finite  element  and  strength  grids  for  concrete  specimens. 
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Table  1.  Description  of  Specimens 


Specimen  Designation 
(L-L-(Xm/ds)] 

L 

(in.) 

d. 

(in.) 

ms 

d, 

(in.) 

(in.) 

3-3-24 

3 

3 

1 

4 

0.25 

6 

3-3-6 

3 

3 

1 

1 

1.00 

6 

12-12-6 

12 

3 

4 

4 

1.00 

6 

The  strength  distribution  functions  were  evaluated  by  generating 
20  realizations  of  X0(x,  y)  over  a  5  x  5-foot  area.  Each  realization  was  equ;va- 
lent  to  casting  a  5  x  5-foot  concrete  slab,  which  could  be  subdivided  into  N 
specimens.  The  ranges  of  x  and  y  provided  an  N  of  25  for  the  specimens  with 
the  largest  size  (12  x  12  inches).  Two  sets  of  25  specimens  of  3x  3  inches  wit'> 
different  strength  grids  were  then  selected  from  the  same  regions  of  15x15 
inches  for  each  realization.  This  was  done  to  determine  the  influence  of  fine¬ 
ness  of  the  strength  grid  when  all  other  -ariables  are  held  invariant. 

Results  and  Discussion 


Figure  12  shows  the  relations  among  the  cumulative  strength  distribu¬ 
tions  determined  from  one  of  the  realizations.  Typically,  the  cumulative 
distributions  were  almost  identical  for  the  two  sets  of  3  x  3-inch  specimens 
that  differed  only  in  the  degree  of  fineness  of  the  strength  grid.  For  each  of 
the  20  realizations,  the  distribution  function  for  the  larger  specimens  (12  x 
12-inch)  fell  to  the  left  of  that  of  the  smaller  specimens.  Table  2  gives  the 
means  and  standard  deviations  obtained  from  each  realization  as  well  as  their 
overall  values  for  all  20  realizations.  The  distribution  functions  of  the  means 
are  shown  in  Figure  13.  This  figure  shows  that  the  fluctuation  of  the  means 
was  smaller  for  the  larger  specimens.  Figure  14  gives  the  overall  strength  dis¬ 
tribution  functions  for  the  3  x  3-inch  and  12x1 2-inch  specimens  as  determined 
from  the  simulation  procedure.  For  reference,  the  idealized  strength  distribu¬ 
tion  function  of  the  reference  specimen  is  also  shown  in  Figure  1 4. 

These  results  suggest  that  the  strength  distribution  function  is  not 
influenced  by  the  fineness  of  the  strength  grid  when  the  ratio  X^/dj  was  taken 
large  enough.  This  result  is  very  significant  because  the  strength  is  an  intrinsic 
property  of  the  material  and  should  not  be  dependent  on  the  strength  grid, 
which  is  an  option  in  the  hand  of  the  analyst.  This  result  could  not  have  been 
obtained  from  the  classical  approach  because  Xm/ds  is  hypothesized  to  be 
always  zero. 


Table  2.  Summary  of  Results 


Means  and  Standard  Deviations  for — 

3-3-24 

3-3-6 

12-12-6 

Specimen 

Specimen 

Specimen 

S25 

°25 

S25 

°25 

^25 

°2S 

(psi) 

(psi) 

(psi) 

(psi) 

tp>-) 

(psi) 

8,036 

1.450 

8.063 

1.434 

5262 

1.019 

8,616 

1,118 

8.637 

1,107 

5273 

1.096 

8,707 

1.346 

8.725 

1.339 

5.806 

1.020 

8,169 

1.394 

8.182 

1.382 

5.678 

1.338 

8,467 

1.020 

8.493 

1.017 

6274 

1.084 

8.414 

1.172 

8.424 

1.171 
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1.409 
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1.665 

8.534 

1.652 

6.403 

c?° 

8.155 
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8.176 

1.064 

6.582* 

899 

8.221 
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1.313 
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8.147 

1.674 

8.169 
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8.753 

846* 

8.775 

846* 

5.674 

1.080 
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1,632 

8,716 

1.623 

6.105 

1.147 

8.669 

1.375 
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1.366 

6.086 

759* 
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1270 

5.366* 

1.693* 
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1.510 

9  411* 
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7.7*6 

1.095 
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9.006 
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1243 
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8282 

1.663 
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1236 
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s500 
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1.614 
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*  Maximum  value 

*  Minimum  value 
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Figure  14.  Cumulative  strength  distributions  from  twenty  realizations, 


The  results  presented  above  followed  the  general  trend  obtained 
from  laboratory  observation  in  that  the  mean  value  of  the  strength  as  well 
as  the  standard  deviation  decreases  as  the  size  of  specimen  increases  (Figure 
14).  The  change,  in  percent,  of  the  mean  strength  with  size  (Figure  14)  is  in 
agreement  with  those  given  in  Reference  15.  However,  one  may  argue  that 
the  values  of  the  parameters  chosen  for  the  numerical  example  cause  such 
agreement.  This  is  true  and,  perhaps,  one  of  the  important  features  of  the 
proposed  approach.  The  ability  to  vary  the  input  parameters  and  compare 
the  results  with  those  obtained  experimentally  may  eliminate  conducting 
new  and  expensive  experiments  to  determine  the  spatial  strength  distribu¬ 
tion  in  the  manner  previously  explained. 

CONCLUSION 

A  new  interpretation  of  the  statistical  size  effect  is  proposed  based 
on  consideration  of  the  spatial  strength  variation  as  a  random  process.  One 
limiting  care  of  the  proposed  approach  is  the  classical  approach  which  assumes 
the  strength  to  be  identically  and  independently  distributed.  The  other  limit¬ 
ing  case  presents  the  engineering  approximation  of  uniform  strength.  The  new 
approach  provides  a  powerful  tool  for  predicting  the  effect  of  size  on  concrete 
strength  in  a  manner  th?<  is  consiste.-f  with  laboratory  observations. 
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LIST  OF  SYMBOLS 


A  Material  constant 

A#  Area  of  an  element 

Am  Maximum  aggregate  size 

a  Argument  of  the  distribution  function 

c  Thickness 

d.  d\  tie  Dimension  of  an  element 

ds  Distance  between  intersection  points  of  a 

strength  grid 

dv  An  element  of  volume  v 

FAla)  Distribution  function  of  the  minimum  value 
of  the  process  X0(x.y)  for  an  element  of  an 
area  Ae 

FA(a)  Distribution  function  of  the  minimum  value 
of  the  process  G(x,  y)  for  an  element  of  an 
area  Ae 

Fd(a)  Distribution  function  of  the  minimum  value 
of  the  process  X0(x)  for  an  element  of 
length  d 

Fd(a)  Distribution  function  of  the  minimum  value 
of  the  process  G(x)  for  an  element  o'  length  d 


Fd(a)  Distribute,..  function  Fd(a)  determined  from 
laboratory  experiment 

Fx(a)  Distribution  function  of  a  random  variable  X 
F„  (a)  Distribution  function  of  a  random  variable  x 


Rc  A  reference  strength  j 

Sn  Mean  strength  of  a  sample  of 

| 

S'20  Mean  of  20  values  of  Sn  for  i 
S20  =  S500>  1 


F 

9 

h 


Homogeneous  Gaussian  process  with  mean 
zero,  unit  standard  deviation,  and  autocor¬ 
relation  function  r 


V  Volume  ^ 

3> 

> 

X  A  random  process  represent! 

of  the  strength  from  the  mai 

XQ  A  random  process  represent! 

the  strength 


Beam  height 


5 

! 

F 

Vr 

p 

is 

Material  constant 
Wave  number 
Maximum  wave  number 
Length  of  a  specimen 
Large  positive  integer 
Positive  integer 
The  ratio  L/de 
The  ratio  d^d. 

Positive  integers 
Positive  integer 

Number  of  intersection  points  of  a  strength 
grid  at  which  the  value  of  XQ(x,  y)  is  evalu¬ 
ated  within  an  element,  e 

Autocorrelation  function 

Autocorrelation  function  r  ew.uated  from 
laboratory  experiment 

A  random  process  describing  the  statistical 
variation  of  strength 

A  reference  st'ength 

Mean  strength  o*  a  sample  of  n  specimens 

Mean  of  20  values  cf  S.,  for  n  =  25  (that  is. 
^20  *  ^500* 

Volume 

A  random  process  representing  the  deviation 
of  the  strength  from  the  mean  Xm 

A  random  process  representing  a  measure  of 
the  strength 


Xm  The  mean  value  of  X0 

Xm  Mean  value  of  Xm  determined  from  laboratory 
experiments 

Xv  Random  variable  representing  the  strength  of 
a  piece  of  material  of  volume  v 

x.  y,  2  Space  coordinate 

A  Dimension  of  a  reference  specimen  or  gage 

length 

5  Dirac  delta  function 

0  A  measure  of  the  stress  path  (=  arctan  o2/oi ) 

Xm  Wavelength  associated  with  km 

£  Distance  between  two  material  points  along 

a  coordinate  axis 

°1  •  °2  Principal  stress  (positive  for  compression) 

0o  Standard  deviation  of  population 

on  Standard  deviation  of  specimen  strength 

for  a  sample  of  n  specimens 

a’20  Standard  deviation  for  20  values  of 

S„  (n  =  25) 
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